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The origin of quark and lepton masses is one of the outstanding problems of physics. 
As the experimental data becomes more and more accurate, testing theories of 
fcrmion masses requires greater care. In this talk we discuss a theoretical frame- 
work for testing those theories with a high energy desert. It is only with precision 
tests that we can hope to narrow the set of viable, beyond the standard model 
theories. 

1 Theoretical framework : top - down approach 

We consider theories in four dimensions with a high energy desert up to a 
physical threshold M* with M pl > M* > M G , the GUT scale. Above M* we 
assume a fundamental theory (such as string theory) from which the effective 
low energy theory is derived with boundary conditions at M*. Below M* we 
have two possibilities - 

• M* 3> Mq and we have a SUSY GUT below. In this case, the relevant 
boundary conditions for low energy phenomenology are defined at Mq. 

• Or we have the minimal SUSY spectrum below M* with relevant bound- 
ary conditions there. 

The latter possibility is the default when talking of string unification. I focus 
on the former possibility in the rest of this talk. 

1.1 SUSY GUT with boundary conditions at Mq 

In this case we have three types of boundary conditions specifying the pa- 
rameters for gauge couplings, Yukawa couplings and soft SUSY breaking at 
M G . 

• Gauge couplings ojj = cxq for i = 1, • • • , 3. 

This is a tree level relation in a GUT. A good fit to the low energy data 
already requires one loop threshold corrections; hence these corrections, 
predominantly from the GUT breaking and Higgs sectors of the theory, 
are observable and thus constrain GUTs. 
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• Yukawa couplings 

Schematically we assume 
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where e, e'/e are small parameters. One generally (with or without 
GUTs) assumes some family symmetry; whether it is a permutation, 
U(l), or non-abelian SU(2), SU(3), or discrete non-abelian symmetry 
D(3). Family symmetries are needed for several reasons: 

— they reduce the number of fundamental parameters; 

— they can explain the hierarchy of fermion masses; 

— they can suppress flavor violation, such as \i — > ej; 

— they are necessary to suppress dimension 5 baryon number violating 
operators responsible for nucleon decay, i.e. they explain the small 
Yukawa-like couplings appearing in these dimension five operators. 
Or they can be used (in theories without GUTs) to forbid, with U(l) 
or baryon parity, dimension 5 baryon number violating operators all 
together; 

— and they can eliminate any other undesirable effects of Planck scale 
physics. 

In general, there are three Yukawa matrices given by 



with (k = u, d, e). The parameters A, B, B, C, D, D are arbitrary, 
order one complex coefficients; fit to the data. 

• SUSY breaking 

For example, in the CMSSM we have too, Ao, tan/3, sign{fi). 

In other SUSY breaking scenaria (such as, GMSB, AMSB, gMSB, ... ) 
we have a different set of SUSY breaking parameters. The main point 
is that simple mechanisms for SUSY breaking have a small number of 
arbitrary parameters; hence they can be predictive. 
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2 Precision Electroweak Data : including fermion masses and 
mixing angles 

The precision electroweak data for testing theories beyond the standard model 
includes, in addition to the electroweak data used for testing the standard 
model, fermion masses and mixing angles. 

• SM (CMSSM) Precision Data 



M z = 91.188 (0.091) 
Gf, ■ 10 5 = 1.1664 (0.0012) 
a s (M z ) = 0.1181 (0.002) 
B(b -» 87) • 10 4 = 2.96 (0.; 

• BSM Precision Data 



M w = 80.419 (0.080) 
a- E x M = 137.04 (0.14) 
Pnew ■ 10 3 = -0.20 (1.1) 



M t = 174.3 (5.1) m b {m b ) = 4.20 (0.20) 

M b -M c = 3.4 (0.2) m c (m c ) = 1.25 (0.10) 

m s (2GeV) = 110 (35) m d /m s = 0.050 (0.010) 
Q = 22.7 (0.8) 

M T = 1.777 (0.0018) Mn = 0.10566 (0.00011) 
M e • 10 3 = 0.5110 (0.00051) 

\V US \ = 0.2196 (0.0023) \V cb \ = 0.0402 (0.0019) 
\V ub /V cb \ = 0.090 (0.025) 

e K ■ 10 3 = 2.28 (0.013) B K = 0.85 (0.15) 

I have included only a subset of the SM precision data. This is the subset 
we have used in our analysis a. Instead of the complete set of standard 
model data we have added the thirteen observables associated with charged 
fermion masses and mixing angles. With regards to the BSM precision data, 
M(m) refer to pole (MS running) masses. M b — M c has smaller theoretical 
uncertainties than m c (m c ). This is because one can use HQEFT to relate the 
difference to Mb — Mr, with renormalon effects cancelling. Hence we have 
used the former rather than the latter. Q is the Kaplan- Manohar- Leutwyler 
ellipse parameter defined by 

2 i 2 

m„ 1 ^ 

vn\ Q 2 m 2 d 

Note, this is the best known relation among the three light quark masses. 
The experimental results X^ p (au) are given above. These are for the most 
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part taken from the PDG 2000. We have however used only the lower half 
of the PDG range for m s ; as suggested by recent lattice results. Finally, in 
those cases where the experimental uncertainty is less than 0.1% we assume 
a numerical error of 0.1% for Ofc. We then define a \ 2 function given by 

N 

2 \ * / v exp v th\2 / 2 

/c=i 

2.1 Global x 2 fits 

We fit the data by minimizing the \ 2 function using the program Minuit. For 
the most part, the parameters we vary are defined at the GUT scale, except 
for /i and tan/3 defined at Mz- We choose boundary conditions at Mq given 
by 

• ax = a.2 = «g, a3 = He (1 + £3) (where £3, parametrizing the one loop 
threshold corrections to gauge coupling unification, is a free parameter); 

• a model dependent parametrization of Yukawa couplings; and 

• soft SUSY breaking parameters which we take to be universal squark 
and slepton masses - m M , universal gaugino masses - M 1 / 2 , universal cubic 
scalar coupling - Aq, non-universal H u , Hd masses, /i and t&n f3(Mz). 

We use two loop renormalization group running from Mq to Mz within 
the MSSM. We include one loop threshold corrections to fermion masses and 
electroweak observables at Mz- Below Mz we use 3 loop QCD + 1 loop QED 
running. Electroweak symmetry breaking is obtained self-consistently using 
the one loop improved Higgs potential, including mf and rai corrections with 
an effective 2 Higgs doublet model defined below M stop 0- Finally the \ 2 
function includes (for three families) 20 low energy observables. 

Our analysis is in the same spirit as the global fits of the standard model 
(or MSSM) to precision electroweak data 13, although we use only a subset flf 
the electroweak data. In the rest of this talk I focus on our recent analysis u. 
including only the third generation quark and lepton masses. 

3 SO(10) Yukawa unification for third generation and 
threshold corrections at Mz 

In the simplest version of SOiq the Yukawa couplings for the third generation 
are all equal at Mq, i.e. A& = At = A T = A„ T = A. The low energy values of 
the top, bottom and r couplings depend solely on the value of Mq, cxq and 
A. The first two are already fixed by gauge coupling unification, hence we 
need only one low energy observable to fix A. We use the ratio mb/m T . We 
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then use m T — \ T -^cosf3 to fix the value of 0; v is fixed by electroweak data. 
Finally we obtain one prediction for nit = At -^g sin/3 ~ 170 ± 20 GeV i. 

3.1 Weak scale threshold corrections - large tan[3 

This beautiful result is however spoiled. by potentially large weak scale thresh- 
old corrections proportional to tan [3 U. By far the largest corrections come 
from SUSY mass insertions given by Srrib ~ 8rrv\ + 8m£ where 5mf oc A "" 8 *°"^ 

b 

and 5m b oc aAil^l§._ These corrections can be as large as ~ 50 %; however 

t 

good fits require Smt < —2%. 

Note that when 5m b > (fj, > 0, in our convention), then Sm b < 0. 
Recent data favors fi > 0. For example, the process b — > sj is obtained by 
similar graphs with a photon insertion. In this case the chargino term typically 
dominates and has opposite sign to the SM and charged Higgs contributions 
for [i > 0, thus reducing the branching ratio. This is a good thing, since the 
standard model plus charged Higgs contribution is typically larger than the 
experimental value. Thus /i < is problematic. Also the sign of the SUSY 
contribution to the anomalous magnetic moment oi the muon is correlated 
with the sign of m,. Recent BNL results favor fi > 0. 

Pierce et al. B and Baer et al. □ find Yukawa unification only when Smf < 
and hence poor fits for b — * sj. On the other hand we find good fits for 
Sm 9 b > H. What is going on: 

In the past, we have neglected the small "wavefunction renormalisation" 
corrections to Snib- These corrections, 5m b ° 9 , are not proportional to tan/3 
nor to At or m~ g . They are only logarithmically dependent on SUSY particle 
masses. They are however positive and about 6%. They thus have to be 
cancelled in order to find good fits to the bottom mass. Note, the corrections 
to the top and r mass are jSmall. In the present analysis □ we include all one 
loop threshold corrections u to quark, charged lepton masses and electroweak 
observables. 

We find Yukawa unification possible only in a narrow region of soft SUSY 
breaking parameter space. 

In this region of parameter space with large tan (3 we need large mie, Aq 
and smaller /i, Mi/ 2 . A large negative value of Aq serves two purposes: 

• it enhances the chargino contribution and 

• it induces large splitting between t\ — t% ■ 

Both effects go towards enhancing Sm^ compared to 8vnP b + Sm b ° 9 . In ad- 
dition large values of m\§ enable the largest possible splitting of ti t % while 
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Figure 1. \ 2 contours for m w = 1500 GeV (Left) and m u = 2000 GeV (Right). The 
shaded region is excluded by the chargino mass limit m + > 103 GeV. 



maintaining a positive stop mass squared. 

This range of parameters is however consistent with electroweak symme- 
try breaking only if the H u , masses are split. We have considered two 
possible ways of splitting the Higgs masses we refer to as Dx term splitting 
and universal scalar masses. 

Dx splitting 

In this case we have rn^ H H ^ — m\ a =p 2 Dx, m (Q « g) = m i6 + Dx and 
TO (d l) ~ TO i6 — 3I5x- Unfortunately, since Dx > we find to| << to| which 
is bad. 

Universal scalar masses 

In the other case we have rn 2 H H s = mf Q =p 2 Ato^, and m 2 Q g ) = 
m (<2 l) = m ?6- W e now find to? << to? which is just what is needed, i.e. 
good. 

In figure 1 we give contours of constant x 2 as a function of [i, M 1 / 2 
with to i6 = 1500 and 2000 GeV. All other parameters have been varied in 
order to minimize \ 2 ■ I n figure 2 we show TO(,(to;,) and <5to{, contours as a 
function of /x, M1/2 with mi6 = 2000 GeV. By comparing the two figures, it 
is clear that x 2 is smallest when Srrib < —2%. The best fits are obtained with 
Aq ~ —1.9 toi6 and toio ~ 1.4 toi6- 

In a more complete theory including all three generations, we would be 
able to calculate the anomalous magnetic moment of the muon a^ EW and 
the proton lifetime. Here we showjesults obtained in a particular 3 family 
model (see Blazek, S.R. and Tobe El). We find a^ EW is less than the recent 
BNL measurement a^ EW = 43 ± 16 x lO" 10 . 1 In fact EW decreases with 
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Figure 2. Contours of constant m;,(m(,)[GeV] (Left) and Smi, in % (Right) for mi6 = 2000 
GeV. 



increasing m\§. With a minimum value of miQ ~ 1200 GeV, necessary for 
reasonable fits, we find a^ EW < 16 x 10~ 10 . 

The proton lifetime, on the other hand, increases as mi6 increases. It can 
be consistent with the most recent Super K preliminary results LLj, but only 
if we allow for the three quark matrix element f3 between a proton state and 
the vacuum to be of order 0.5 x/3; attice IB. 



4 Conclusion 



In conclusion, theories beyond the standard model are strongly constrained 
by precision electroweak data, which now includes fermion masses and mixing 
angles. As an example, we have discussed SO(10) Yukawa unification. We 
find that good fits constrain the soft SUSY breaking parameters to a narrow 
region given by mi6 > 1200 GeV, A® ~ —1.9 mi6, and mio ~ 1-4 m\Q. This 
then implies significant predictions for Higgs and SUSY particle masses. See 
the talk of R. Dermfsek in these proceedings for further details. 
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Figure 3: and the proton lifetime plotted as a function of fi, M1/2 for mi6 = 

2000 GeV. 
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